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Abstract 

In the present paper, an exact analytic solution for the optimal unambiguous state 
discrimination(OPUSD) problem involving an arbitrary number of pure linearly indepen- 
dent quantum states with real and complex inner product is presented. Using semidefinite 
programming and Karush-Kuhn- Tucker convex optimization method, we derive an ana- 
lytical formula which shows the relation between optimal solution of unambiguous state 
discrimination problem and an arbitrary number of pure linearly independent quantum 
states. 
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Many applications in quantum communication and quantum cryptography are based on trans- 
mitting quantum systems that, with given prior probabihties, are prepared in one from a set of 
known mutually nonorthogonal statesflj. A fundamental aspect of quantum information theory 
is that nonorthogonal quantum states cannot be perfectly distinguished. Therefore, a central 
problem in quantum mechanics is to design measurements optimized to distinguish between a 
collection of nonorthogonal quantum states. The topic of quantum state discrimination was 
firmly established in the 1970s by the pioneering work of Helstrom [2], who considered a min- 
imum error discrimination of two known quantum states. In this case, the state identification 
is probabilistic. Another possible discrimination strategy is the so-called unambiguous state 
discrimination (USD) where the states are successfully identified with nonunit probability, but 
without error. 

USD was originally formulated and analyzed by Ivanovic, Dieks, and Peres [31 IH E] iii 
1987. The solution for unambiguous discrimination of two known pure states appearing with 
arbitrary prior probabilities was obtained by Jaeger and Shimony [6]. Although the two- 
state problem is well developed, the problem of unambiguous discrimination between multiple 
quantum states has received considerably less attention. The problem of discrimination among 
three nonorthogonal states was first considered by Peres and Terno [5]. They developed a 
geometric approach and applied it numerically on several examples. A different method was 
considered by Duan and Guo ^ and Sun et al.[fij . Chefles [Qj showed that a necessary 
and sufficient condition for the existence of unambiguous measurements for distinguishing 
between N quantum states is that the states are linearly independent. He also proposed a 
simple suboptimal measurement for unambiguous discrimination for which the probability 
of an inconclusive result is the same regardless of the state of the system. Equivalently, the 
measurement yields an equal probability of correct detection of each one of the ensemble states. 
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Over the past years, semidefinite programming (SDP) has been recognized as a valuable 
numerical tool for control system analysis and design. In SDP, one minimizes a linear func- 
tion subject to the constraint that an affine combination of symmetric matrices is positive 
semidefinite. SDP has been studied (under various names) as far back as the 1940s. Subse- 
quent research in semidefinite programming during the 1990s was driven by applications in 
combinatorial optimization [10], communications and signal processing [IH [121 [13] and other 
areas of engineering [13]. Although semidefinite programming is designed to be applied in 
numerical methods, it can be used for analytic computations, too. In the context of quantum 
computation and quantum information, Barnum, Saks, and Szegedy have reformulated quan- 
tum query complexity in terms of a semidefinite program [15]. The problem of finding the 
optimal measurement to distinguish between a set of quantum states was first formulated as 
a semidefinite program in 1974 by Holevo, who gave optimality conditions equivalent to the 
complementary slackness conditions [2]. Recently, Eldar, Megretski, and Verghese showed that 
the optimal measurements can be found efficiently by solving the dual followed by the use of 
linear programming [16]. Also in Ref. [17], SDP has been used to show that the standard algo- 
rithm implements the optimal set of measurements. All of the above mentioned applications 
indicate that the method of SDP is very useful. The reason why the area has shown relatively 
slow progress until recently within the rapidly evolving field of quantum information is that it 
poses quite formidable mathematical challenges. Except for a handful of very special cases, no 
general exact solution has been available involving more than two arbitrary states and mostly 
numerical algorithms are proposed for finding optimal measurements for quantum-state dis- 
crimination, where the theory of the semidefinite programming provides a simple check of the 
optimality of the numerically obtained results. In Ref. [H] we obtained the feasible region in 
terms of the inner product of the reciprocal states and the feasible region in terms of the inner 
product of the states which enables us to solve the problem without using reciprocal states. 
Moreover, for the real inner product of states, we obtained an exact analytic solution for 
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OPUSD problem involving an arbitrary number of pure linearly independent quantum states 
by using KKT convex optimization method. In this paper, an exact analytic solution for the 
optimal unambiguous state discrimination(OPUSD) problem involving an arbitrary number of 
pure linearly independent quantum states with real and complex inner product of states,using 
the Karush- Kuhn- Tucker convex optimization method, is given. 

The organization of the paper is as follows. First, the definition of the unambiguous quan- 
tum state discrimination is given. Then, Semidefinite programming, Karush-Kuhn- Tucker 
(KKT) theorem and SDP formulation of unambiguous discrimination is studied. Finally for 
the real and complex inner product of reciprocal states, an exact analytic solution for OPUSD 
problem involving an arbitrary number of pure linearly independent quantum states is pre- 
sented by using KKT convex optimization method. The paper ends with a brief conclusion. 

2 Unambiguous quantum state discrimination 

In quantum theory, measurements are represented by positive operator valued measures (POVMs). 
A measurement is described by a collection of measurement operators. These operators are 
acting on the state space of the system being measured. The index k refers to the measurement 
outcomes that may occur in the experiment. In quantum information theory the measurement 
operators are often called Kraus operators [19]. If we define the operator 



the probability of obtaining the outcome k for a given state pi is given by p{k\i) = Tr{IlkPi). 
Thus, the set of operators 11^ is sufficient to determine the measurement statistics. 




2.1 Definition of the POVM 

A set of operators {11^} is named a positive operator valued measure if and only if the following 
two conditions are met: (1) each operator 11^ is positive positive 4^ (iplUklip) > 0, V|'0) and 
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(2) the completeness relation is satisfied, i.e., 



(2.2) 



k 



The elements of {H^} are called effects or POVM elements. On its own, a given POVM {11^} 
is enough to give complete knowledge about the probabilities of all possible outcomes; mea- 
surement statistics is the only item of interest. Consider a set of known states pi,i — 1, N 
with their prior probabilities rji. We are looking for a measurement that either identifies a state 
uniquely (conclusive result) or fails to identify it (inconclusive result). The goal is to minimize 
the probability of inconclusive results. The measurements involved are typically generalized 
measurements. A measurement described by a POVM {Ilk}^^^ is called an unambiguous state 
discrimination measurement (USDM) on the set of states {pi}^i if and only if the following 
conditions are satisfied. (1) The POVM contains the elements {^k}k=o where N is the number 
of different signals in the set of states. The element IIo is related to an inconclusive result, 
while the other elements correspond to an identification of one of the states pi,i — 1,...,N. 
(2) No states are wrongly identified, that is, Tr{piUk) ^ ^ k,k ^ 1, ...,N. Each USD 
measurement gives rise to a failure probability, that is, the rate of inconclusive result. This 
can be calculated as 



A measurement described by a POVM {11^^ } is called an optimal unambiguous state discrim- 
ination measurement (OUSDM) on a set of states {pi} with the corresponding prior prob- 
abihties {rji} if and only if the following conditions are satisfied. (1) The POVM {H^^*} is 
a USD measurement on {pi}. (2) The probability of inconclusive result is minimal, that is, 
(5({n^^*}) = min Q{{Ilk}), where the minimum is taken over all USDM. Unambiguous state 




(2.3) 



The success probabihty can be calculated as 




(2.4) 
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discrimination is an error-free discrimination. This implies a strong constraint on the mea- 
surement. Suppose that a quantum system is prepared in a pure quantum state drawn from 
a collection of given states {iT/^j)},! < i < N in d-dimensional complex Hilbert space "H with 
d > N . These states span a subspace U of "H. In order to detect the state of the system, 
a measurement is constructed comprising + 1 measurement operators {IljjO < i < N}. 
Given that the state of the system is I'f/'j), the probability of obtaining outcome k is {il)i\U.k\'4^i) ■ 
Therefore, in order to ensure that each state is either correctly detected or an inconclusive 
result is obtained, we must have 

(V'iinfclV^i) = pi5ij, l<i,k<N (2.5) 

for some < Pi < 1. Since Uq = Id — J2iLi^iy have {i/jilUoli/ji) = 1 — pj. So a system 
with given state l^/^j) , the state of the system is correctly detected with probability pi and an 
inconclusive result is obtained with probability 1 — pi. It was shown in Ref. [20] that Eq.f l2.5p 
is satisfied if and only if the vectors {i/ji) are linearly independent, or equivalently, dim U = N . 
Therefore, we will take this assumption throughout the paper. In this case, we may choose 

m 

Ii,=Pi\^i){iPi\, l<i<N, (2.6) 

where the vectors I'i/'j) G U are the reciprocal states associated with the states , i.e., there 
are unique vectors in U such that 

= 5ij, l<i,k<N. (2.7) 

With $ and $ we denote the matrices such that their columns are \ilJi) and respectively. 
Then, one can show that $ is 

$ = $($*$)-!. (2.8) 

Since the vectors = 1, are linearly independent, $*$ is always invertible. Alterna- 

tively, 

$ = ($$*)*$, (2.9) 
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so that 

|^,) = ($r)*|^,), (2.10) 

where (...)"'■ denotes the Moore-Penrose pseudoinverse [22l [23] . The inverse is taken on the 
subspace spanned by the columns of the matrix. If the state {ipi) is prepared with prior 
probabihty r/j, then the total probability of correctly detecting the state is 

N N 

p=j2vi{H^i\A) = J2'^iPi (2.11) 

i=l i=l 

and the probability of the inconclusive result is given by 

N N 

Q = l-P = J2vi{i^i\Wi) = 1 - (2.12) 

1=1 i=l 

In general, an optimal measurement for a given strategy depends on the quantum states and 
the prior probabilities of their appearance. In the unambiguous discrimination for a given 
strategy and a given ensemble of states, the goal is to find a measurement which minimizes the 
inconclusive result. In fact, it is known that USD (of both pure and mixed states) is a convex 
optimization problem. Mathematically, this means that the quantity which is to be optimized 
as well as the constraints on the unknowns, are convex functions. Practically, this implies 
that the optimal solution can be computed in an extremely efficient way. This is therefore 
a very useful tool. Nevertheless, our aim is to understand the structure of USD in order to 
relate it with neat and relevant quantities and to find feasible region for numerical and analytic 
solutions. 

3 Semidefinite programming 

A SDP problem requires minimizing a linear function subject to a linear matrix inequality 
(LMI) constraint 

minimize p = (Fx, such that F(x) > 0, (3.13) 
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where is a given vector, x — {xi, ...,Xn), and F{x) — Fo + X^^XjFj, for some fixed Hermitian 
matrices Fj. The inequality sign in F{x) > means that F{x) is positive semidcfinitc. This 
problem is called the primal problem. Vectors x whose components are the variables of the 
problem and satisfy the constraint F{x) > are called primal feasible points, and if they 
satisfy F{x) > they are called strictly feasible points. The minimal objective value c^x is by 
convention denoted by P* and is called the primal optimal value. Due to the convexity of set 
of feasible points, SDP has a nice duality structure with the associated dual program being 

maximize -Tr[FoZ], Z >, Tr[FiZ] ^ Ci. (3.14) 

Here the variable is the real symmetric (or Hermitian) matrix Z, and the data c, Fj are the 
same as in the primal problem. Correspondingly, matrix Z satisfying the constraints are called 
dual feasible (or strictly dual feasible if Z >). The maximal objective value of — Tr[FoZ], i.e., 
the dual optimal value, is denoted by d*. The objective value of a primal (dual) feasible point 
is an upper (lower) bound on P*{d*). The main reason why one is interested in the dual 
problem is that one can prove that d* < P*, and under relatively mild assumptions, we can 
have P* = d*. If the equality holds, one can prove the following optimality condition on x. A 
primal feasible x and a dual feasible Z are optimal which is denoted by x and Z if and only if 

F{x)Z = ZF{x). (3.15) 

This latter condition is called the complementary slackness condition. In one way or another, 
numerical methods for solving SDP problems always exploit the inequality d < d* < P* < P, 
where d and P are the objective values for any dual feasible point and primal feasible point, 
respectively. The difference 

P*-d* = Jx + Tr[FoZ] = Tr[F^Z] > (3.16) 

is called the duality gap. If the equality holds d* = P*, i.e., the optimal duality gap is zero, 
then we say that strong duality holds. 
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3.1 Karush-Kuhn-Tucker (KKT) theorem 

Assuming that functions gi, hi are differentiable and that strong duahty holds, there exists 
vectors ^ E R'^ and y e such that the gradient of dual Lagrangian L{x*.^*,y*) — f{x*) + 
+ '^iViOii^*) ^^^^ ^ vanishes at x*: 

hi{x*) = {primal feasible), 

gi{x*) < {primal feasible), 
y* > {dual feasible), ylgi{x*) = 0, 
V f(x*) + 5] C V hi{x*) + J] y* V 9i(x*) = 0. (3.17) 

i i 

Then x* and {Ci,y*) are primal and dual optimal with zero duality gap. In summary, for any 
convex optimization problem with differentiable objective and constraint functions, the points 
which satisfy the KKT conditions are primal and dual optimal, and have zero duality gap. 
Necessary KKT conditions satisfied by any primal and dual optimal pair and for convex prob- 
lems, KKT conditions are also sufficient. If a convex optimization problem with differentiable 
objective and constraint functions satisfies Slaters condition, then the KKT conditions provide 
necessary and sufficient conditions for optimality: Slaters condition implies that the optimal 
duality gap is zero and the dual optimum is attained, so x is optimal if and only if there are 
{^ijVi) such that they, together with x, satisfy the KKT conditions. 

3.2 Slaters condition 

Suppose X* solves 

minimize f{x)gi{x) >bi, i = l,...,m (3.18) 
and the feasible set is nonempty. Then there is a nonnegative vector ^ such that for all x 

L{x,0 = f{x)+e[b^g{x)] < f{x*) = L{x*,0- (3.19) 
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In addition, if f {...), gi{---), = 1, are continuously differentiate, then 

f^-?^ = 0. (3.20) 

o{xj) d{x) 

In the spatial case the vector x is a solution of the linear program 

minimize c^x, such that Ax = bx > 0, (3-21) 

if and only if there exist vectors ^ & and y G for which the following conditions hold 

for {x,^,y) = {x\C,y*)- 

A^^ + y = c,Ax = b,Xi>0; y > 0; 

Xiyi = 0, i = l, ...,m. (3.22) 

A solution {x*,^*,y*) is called strictly complementary, if x* + y* > , i.e., if there exists no 
index i G 1, ...m such that x* = y*. 



4 SDP Formulation of unambiguous discrimination 

Eldar, Megretski, and Verghese in Ref. [21] have showed that the unambiguous discrimination 
problem can be reduced to SDP method and the KKT conditions can be defined as 

N 



F{P) = J2p'^' + Fo>0, Tr{UiX) = Zi + r]i 
1=1 

N 



> 0, 1 < z < iV, AF{p) = 0, X{Ia - J2P^^^) = 0' 

i=l 

N 

ZiP^ = 0, 1 < i < iV, : ^PiFi + Fo > 0, 



i=l 



3X,;z : X > 0,2 > (4.23) 
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where is the x matrix and Y is the N x d matrix. 
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(4.24) 



5 Analytical solution of unambiguous state discrimina- 
tion problem 

When the set of states {ipi}fLi and the prior probabihties {rii}fLi are given, the POVM elements 
for the measurement (Ilj = J^iLiPili^dii^ili = / — J2iLi Yli) depend on only the variables 
Pi and these variables lie in the feasible region. Feasible region that gives the domain of 
acceptable values of pi is determined by the following equation [251 ISSl [271 128] : 

1 - X] DiPi + DijPiPj - DijkPiPjPk + •■• + (-1)^ Y PiiPi2-PiNDni2...ir, = 0, 

i i<j i<j<k ii<i2...<ijv 

(5.25) 
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where the set of {Aii2...«]v} the subdeterminants (minor) of matrix D defined by 

\ 



D 



^ fill 5.12 ■ ■ ■ O^lN 
fi21 0,22 ■ ■ ■ 0'2N 



(5.26) 



\ om aN2 ■ ■ ■ j 
with aij = 

In this section, for the real and complex inner product of reciprocal states, an exact ana- 
lytic solution for OPUSD problem involving an arbitrary number of pure linearly independent 
quantum states is presented by using KKT convex optimization method. The KKT conditions 
for unambiguous discrimination of N linearly independent states are given by 

N 



i=l 



N 



N 



[I - 5^Pd^^)(V^.|)x = x(/ - J2p^\^^)(^^\) = 0' 



1=1 



i=l 



ZiVi = 0, 



Tr{X\^i){tPi\) = Zi + ri„r], > 0, ^ = 1, N 



(5.27) 



Then using KKT conditions one can show that 



1 - piaii -J92«12 
-^1^21 1 - ^2^22 



-Pnclin 

—PnO'2N 



\ 



( 



X\\ X\2 
^21 ^22 



XiN 
X2N 



(5.28) 



\ —Pl^Nl —P2ON2 • • • 1 — PN^NN ) \ ^Nl XN2 ■ ■ ■ XmN J 

If Pi > 0, (i = 1, A^) usin Eg. (15.27]) we can conclude that X is the rank one matrix and we 
have 

x,^ = e^'^'e^^^V^^V^ (5.29) 
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Consequently, the Eq. f l5.28p can be written as 



1 -piau 


-P2ai2 


-PNaiN ^ 








1 - ^2^22 ■ ■ 


—PNl2N 


X 




-Piam 


-P2aN2 ■ ■ 


■ 1 - PNaNN ) 




\ XiN J 



Or, equivalently. 



(5.30) 



^ an ai2 



^21 ^22 



aiN 

tt2N 



X 



Xii 
a;22 








P2 



Xl2 



(5.31) 



\ am aN2 ■■■ J \ ■ ■ ■ x^j^ } \ Pm / \ ) 

If the first condition in Eq. fl5.27p is satisfied, using above equation one can show that the 
optimal solution for pi is given by 



Vi 



e'f-^idetiD) 



X det 



02,1 



a2,j-i e"^^y/r]N a2,i+i 



0.2,N 



(5.32) 



\ cln,i ■■■ CLN,i-i e'^'y/r]i ajy^j+i ■■■ aN,N ) 
Using the fact that the matrix whose components are is inverse matrix of the one whose 
components are aij = (ipilipj) , we have 



XllPl 

X12P2 



\ 



( 



0^21 0-22 



aiAT 

027V 



X\2 



XiiPi 



N 

E 



^ij X{j 



\ XinPn ) \ O-Nl ClN2 ■ ■ ■ ClNN / \ XiN J 

Using Eq. (I5.29P we have 
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Then we can rewrite the probabihty pi in a simpler form 



Pi 



-lip, 



N 
i=l 



(5.33) 



where e '"^^ajj is determined such that Pi satisfy the condition < Pi < 1. Then the success 
probabihty with respect to dij components is given by 



P 



det(D) 



X det 







\ 



D 



(5.34) 



The success probabihty with respect to aj, components is as fohows 



N 



N N 



'\/ViVj^ij 



1=1 



(5.35) 



If the first condition of KKT is not satisfied, for a specified i, Pi does not he in the feasible 
region {pi < or pi > 1). In this case, one can omit the jth row and jth column in the square 
N X N matrices and jth row in the row matrices in Eqs. f l5.32p and f l5.34p , and find the 
optimal solutions with pi = 0. After this reduction, if any other pi does not lie in the feasible 
region, the same procedure will be repeated. Noted that a similar result was derived in Ref. 
p9] with another method. The equations (I5.33P and fl5.35p gives an analytical relation between 
the maximum average success probability and the pure linearly independent quantum states 
to be discriminated. In general, finding an exact analytic solution for the OPUSD problem 
involving an arbitrary number of pure linearly independent quantum states is hard, since the 
explicit expressions of the phases e^'^^{j = 1, N) are not given in equations fl5.33p and fl5.35p . 
However using f l5.33p and fl5.35p one can simplify the calculation of the optimal solution in 
special cases and it may also drive some bounds for the average success probability [HI [29] . 
Then, the approximated methods are useful for unambiguous discrimination of N linearly 
independent quantum states. Since we have presented an analytic relation for the feasible 
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region of N linearly independent quantum states, and this region is convex, then one can easily 
obtain the optimal POVM by some well-known numerical methods such as constrained linear 
or nonlinear least-squares, interior points, and simplex and quadratic programming methods 

m- 

6 Conclusion 

In conclusion, for the real and complex inner product of states, we have been able to obtain 
an exact analytic solution for OPUSD problem involving an arbitrary number of pure linearly 
independent quantum states by using KKT convex optimization method. Moreover, Using 
semidefinite programming and Karush-Kuhn- Tucker convex optimization method, we have 
been able to obtain an analytical formula which shows the relation between optimal solution 
of unambiguous discrimination problem and an arbitrary number of pure linearly independent 
quantum states to be identified. Using this analytical formula one can simplify the calculation 
of the optimal solution in special cases and it may also drive some bounds for the average 
success probability. 
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